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Applications of the Dirac equation with an anomalous magnetic mo- 
ment are considered for description of characteristics of electrons, muons 
and quarks. The Dirac equation with four-dimensional scalar and vec- 
tor potentials is reduced to a form suitable for a numerical integration. 
When a certain type of the potential is chosen, solutions can approx- 
imate quark states inside hadrons. In view of complicated behaviour 
of quarks in a confinement domain some generalizations are considered 
such as the Dirac-Gursey-Lee equation, the Dirac equation in a five- 
dimensional Minkowski space, the Dirac equation in a quantum phase 
space. Extended symmetries for the Dirac equation and its generaliza- 
tions are considered, which can be used for investigation of properties 
of solutions of these equations and subsequent applications in particle 
physics. 



1. Introduction 

It is well known the Dirac equation for a relativistic electron is one of the 
main building blocks of the quantum electrodynamics (QED) and the Standard 
Model (SM) [1, 2]. It is used to check the validity of SM as to search for new 
physics beyond SM. Although the Dirac equation was offered to the descrip- 
tion of electron it describes with success muons, tayons and moreover quarks. 
The fact of validity of the Dirac equation (DE) for description of quarks is as- 
tonishing because quarks do not observe in our physical spacetime. DE is used 
in the quantum chromodynamics (QCD) [2] and in phenomenological hadron 
models, for example, in the relativistic model for quasi-independent quarks 
(RMQIQ), which has been applied for the description of hadron properties [3, 
4]. However in a domain of confinement the deviations from the Dirac equa- 
tion and standard spacetime symmetry can take place. From this point of view 
generalizations of DE are welcome. We consider some generalized equations of 
Dirac tj^e, namely, the Dirac-Gursey-Lee equation (DGLE), the Dirac equa- 
tion in a five- dimensional Minkowski space (FDDE), the Dirac equation in a 
quantum phase space (QPSDE). 
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2. The Dirac equation and AMMs of electrons, muons and quarks 

In 1928 Dirac presented the electron equation, which is used successfully 
for description of fermions at present. Let us write DE in the natural system 
of units: c = h = 1 (p* = id\ = {t, x}, 7* = j3a}). 

(7>.-mMx) = 0, ^=(^; 0), (1) 

where o"o = l2- o" are the known Pauli matrices, the metrical tensor i^ij — 
diag{l, —1, —1, —1}. DE is invariant under transformations of Poincare group, 
which are generated by operators of momenta and four-dimensional rotations 

i 1 

Jij — XiPj XjPi + Sij, Sij — 2*^*-'' ^^-^ — 2^'^''^-' ^j^i)' (^) 

If we define the pseudoscalar matrix 75 = i7°7^7^7^, the left and right 
wave functions V^^, iI^r can be written 

One can consider the charge conjugated fields = C{ijj)'^, where C is the 
matrix of the charge conjugation which satisfies the relations: C^jJC^^ = —'~fi, 

= —C. We can construct two Majorana fields Xi and X2 using ip and ip'': 
Xi = (V' + ■0'^) / V^: X2 — {ip — '4>'^) I a/S- For the Majorana fields the following 
relation holds (xi,2) = ±(xi,2)'^- 

It is known that for a particle "a" with charge g, spin S*, and mass m the 
magnetic moment arises with the absolute value = g{\q\j'2m)S ^ where g is 
the gyro-magnetic factor or the g-factor. The g-factor is usually written as 
g — g'o(l + o), o is the reduced anomalous magnetic moment or anomalous 
magnetic moment (AMM) of the particle. For electrons and muons g^ — l and 
DE with AMM have the so-called Pauli term 

{1% -m- = 0, (4) 

where p'i = Pi- qAi, = diAj - djAi. 

In the SM framework AMM for electrons and muons have been evaluated 
with high accuracy. It is interesting that the most stringent test of SM and 
QED is the comparison of experimental and theoretical values of AMM for 
electrons (ce). For experimental determination of Og the frequencies of radi- 
ation are measured when an electron transits between quantum levels in a 
cylindrical Penning trap [5]. The QED corrections entered into a value Ug are 
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evaluated up to 10th order (the total contribution of terms of 10th order are 
only estimated now) ) [6] . If the validity of SM is assumed, the most accurate 
value of the fine structure constant can be obtained at present (0.37 ppb). 

g/2{exp) = 1.00115965218073(28), a'^ = 137.035999084(51) (5) 

The analogous calculations have been made for AMM of a muon [7, 8] and have 
been compared with the result of E821 experiment at Brookhaven National 
Laboratory (BNL) [9]: a^,{exp) = 11659208.9(6.3) x 10"^°. When the SM 
prediction a^(S'M) = 11659177.3(4.8) x 10^^° have been compared to the BNL 
measurement result the difference is obtained a^{exp) —af^{SM) = 31.6(7.9) x 
10"-'^°, which corresponds to the Aa discrepancy! To establish physics beyond 
SM new 9^ — 2 experiments are planned [10]. The role of electromagnetic 
AMMs of quarks should be make clear in the future. 

3. The Dirac equation and the relativistic quasi-independent quark 
model 

The main statement of the relativistic model for quasi-independent quarks 
is that hadron can be described as a system of independent constituents (or 
quasi-independent ones with weak residual interactions), which move in some 
mean self-consistent field. In the framework of RMQIQ a wave function ipi 
for any valent i-quark is a sohition of DE with the mean field static potential 
U{ri). Energies or mass terms Ei{nl,ji) represent the energies of constituents 
in a mean field {n\ and j'j are the radial quantum number and the quantum 
number of the angular moment correspondingly for the i— th constituent) and 
can be evaluated for quarks with the help of solutions of the stationary Dirac 
equation (for diquarks and constituent gluons with the help of solutions of the 
Klein-Gordon- Fock equation): 



with Vo(r) = (Tr/2 and Vi{r) — — 2Q;s/3r, where the model parameters a and 
as have meanings of the " string tension" and the strong coupling constant at 
small distances, correspondingly [3, 4]. 

Angular dependence of the single-particle wave functions in a stationary 
state for the spherically symmetric potential can be separated in a well-known 
manner, namely, the solutions of Eq.(6) with the total angular momentum j 
and its projection m can be represented as (the subscript i here and below is 
omitted) 



(6) 




(7) 



3 



where n —r/r. Constituent energies Ei are evaluated by solving the radial 
equations with the model potentials for each constituent. If A; = —oo{j + 1/2), 
where u is connected with an eigenvalue of the space-parity operator, the 
system of the radial Dirac equations reads 

{E - Vo-V,-m)f ^ -^^^—^g - g\ 

r 

{E-rV,-V^+ m)g = -^^^/ + /', (8) 

One can derive the second order equation for the "large" component /(r) 
using Eqs.(8), then making two substitutions E — ^JX-\-■m? and 

0(r) = rf{r)[Vo{r) - T4(r) + m + VA + m2]-V2, 

one comes on to the model radial equation for 0(r) in the following form [3, 
4]: 

0" + Ac/) = [V^ - + 2{mVo + yXT^yi) + M^±ll 



4( VA + m? - + H + m)2 r(VA + m2 - + + m)^ 
2( VA + m? - Vi + H + m)2^ 



(9) 



This equation is suitable for numerical integration and fitting of hadron mass 
spectra. For instance, one can use the phenomenological mass formula for 
g'^-mesons with the quark and antiquark energies, which takes into account 
spin-spin interaction between quark and aniquark, in the following form [3, 4]: 

M(g'g) = Eq + Ei + E2 + A< S1S2 > ^2 (10) 

It is interesting to relate the model of quasi-independent quarks to the 
constituent quark model, where the Zeldovich-Sakharov formula are frequently 
used [11, 12]: 

M(g'g) = mo + mi 7712+ < S1S2 > ^12, (H) 

where mi, m2 are masses of constituent quarks, mo is some additional phe- 
nomenological contribution. Take into account the similarity between formulae 
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(10) and (11) we can relate constituents' energies Ei with constituent quarks' 
masses within uncertainties of models. Some complications arise for exited 
hadron states. In the model of quasi-independent quarks Ei remains an energy 
of i-th exited constituent, while in the constituent quark model an additional 
contribution representing energy of excitation must be taken into account. 

The values of parameters a and ctg entered into potentials Vo(r) and Vi{r) 
have been determined after the comparison between experimental and evalu- 
ated spectra of the vector and pseudoscalar mesons. In this manner it was 
found that the value of a is the same for the light and heavy mesons within 
the systematic errors of the model: a — (0.20 ± 0.01)GeV'^, whereas the ag 
values run from 0.7 to 0.25 for the light and heavy mesons, correspondingly 
[3, 4]. The values of constituents' energies Ei for u—,d—,s—, c— ,6— quarks 
and antiquarks in MeVs are: Eu = 335 ±2, Ed = 339 ±2, Eg = 485 ± 8, 
Ec = 1610 ± 15, Eb = 4952 ± 20 [13, 14]. So RMQIQ may be considered as 
the workable generahzation of the constituent quarks model. In the RMQIQ 
framework the large number of meson resonance masses have been evaluated, 
which do not contradict as the existing experimental data, as the results of 
other phenomenological models [4, 15, 16, 17]. In some variants of the quark 
model anomalous chromo-magnetic moments of quarks have been taken into 
account as well [18, 19, 20]. 

4. The DGL equation and the Dirac equation in a five-dimensional 
Minkowski space 

It is mentioned above, the spacetime symmetry for quarks can differ from 
the Poincare symmetry. Let us consider some extended kinematical symmetries 
and generalized Dirac equations. 

The structure of the algebra of Dirac 7-matrices (including the 75 -matrix 
anticommutating with other 7-matrices) shows that the extension of the Poincare 
group in a four-dimensional Minkowski space to a generalized Poincare group 
in a five- dimensional Minkowski space is possible. By this way the Dirac equa- 
tion is obtained in a five-dimensional Minkowski space [21]: 



where = = {t,x, x^}, 7" = {13, 13a,i^5}, the metrical tensor riah ~ 

diag{l, —1,-1,— 1,-1}. FMSDE is invariant under transformations of Poincare 
group in a five-dimensional Minkowski space P(l,4), which generate by oper- 
ators of momentum components and five- dimensional rotations 



(7X - k)iP{x) = 0, 



(12) 



J 



ah 



-i^ilalb-lbla)- (13) 
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The mass of a particle is dependent on a value of an additional momentum in 
a new direction p^, = + (p^)^, where k is the constant entered in Eq. 

If the symmetry group of a equation was restricted to the homogeneous 
pseudoorthogonal group 0(1,4) (or 0(2,3)), then the equation is depended 
on generators of the pseudoorthogonal rotations Lab and have the form 

l/27VW(x) = A^(x) (14) 

This equation is proposed by Dirac in 1935 [22]. However A has an imaginary 
term due to nonhermitean operators Ljj for a finite dimensional representation. 
The imaginary term should be single out for a correct physical interpretation, 
and the equation (14) was presented by Gursey and Lee in another form [23]: 

{Ai + 7V^afe)^?(a;) = 2mRri{x) (15) 

where A = rriR — 2i, R is a. radius of a de Sitter space. Properties of solu- 
tions of equations of this type and its application for a description of elemen- 
tary particles were studied by Fronsdal, who used the generalized momentum 
Po + dpoL = Pfi where po and L have the forms of the usual generators of trans- 
lations and Lorentz transformations in Minkowski spacetime, and p^ = mg, 
d = fis/mo [24]. Take into account the results obtained in Refs. [22, 23, 24], 
we can apply the following modified Dirac-Gursey-Lee equation for a descrip- 
tion of quark characteristics [25]: 

[7,(p^ + dp'^.Ll + tfisY/2) + 2ifisS,,{U^ + S'^)]tP = rmp. (16) 

To estimate a /is value with the help of a constituent quark mass m and a 
current quark mass mo values, a quark ground state ipo can be used in a 
meson so the contribution from LipQ can be neglected. In this way we obtain 
from Eq.(23) the approximate relation: m ^ mo + 2ifis- So the constant fig 
should be pure imaginary and \fi\s ~ O.lGGeV [25]. 

The problem of quark and lepton masses is the most important problem 
in particle physics now. In the SM framework these masses arise due to the 
Yukawa coupling between fundamental fermions and Higgs field. But, it is 
known the Higgs meson has not been discovered [26] and LHC experiments in 
the near future must make clear the question of its existing. The mass problem 
of quarks and leptons became more complicated after discovering the neutrino 
masses [26]. The neutrino masses are very small, so their origin can differ from 
the origin of other fundamental fermion masses due to the fact that neutrinos 
can be Majorana particles. It is interesting that values of quark and neutrino 
masses can be connected with values of there mixing angles [14]. 
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5. The Dirac equation in a quantum phase space 

It is well known that the Poincare symmetry, which is the spacetime sym- 
metry of the relativistic quantum field theory, originates from the isotropy and 
homogeneity of Minkowski spacetime and is based on observations of macro- 
and micro-phenomena concerning conventional physical bodies and particles. 
However, description of color particles over the whole range of interaction dis- 
tances is possibly needed using extended symmetries and theories with new 
fundamental physical constants (other than the well known ones c and h) in a 
noncommutative spacetime [27]. 

Let us consider the generalized model for a color particle motion, when 
coordinates and momenta are on equal terms and form an eight dimensional 
phase space: h ={h^\h^ = q^'^A = 1,2,3,4, /i = 0,1,2,3, = Tp^',A = 
5, 6, 7,8, j2 = 0, 1, 2, 3}, P = {P^\P^ = p^^A= 1, 2, 3, 4, ^ = 0, 1, 2, 3, P^ = 
aq^, A = 5,6,7,8, n = 0, 1, 2, 3} [28]. The constants r and a have dimensions 
of length and mass square, correspondingly. Their values can be chosen on the 
phenomenological ground or with the help of some functions of the quantum 
constants /i, k and A. We assume the generalized lenght and mass squares 

= h^hA, = P^Pa (17) 

are invariant under the 0(2,6) transformations, where Ha = QAsh^, Qab = 
g^^ = diag{l, -1, -1, -1, 1, -1, -1, -1}. 

Thus the generalized differential mass squared can conserve for strong in- 
teracting color particles: 

dM^ = {dp^f - {dp,f - {dp^f - {dp^f + a\dqof- 

a\dqif - a\dq2f - a\dq^f = {dmf + a\dsf. (18) 

An important point is that the coordinates q^ and the momentum components 
p'^ arc the quantum operators satisfied the generalized Snyder- Yang algebra 
(GSYA) [29, 30, 31, 32]: 

[Fij, Fki] — i{gjkFii — QikFji + guFjk — QjiFik), 

[Fij ,Pk]=i {gjkPi - 9ikPj ) , [Fij ,qk]=i igjkQi - gikQj ) , 

[Fij, I] = 0, \pi, qj] = i{gijl + KFij), (19) 

\pi, I] = i{i/qi - Kpi), [qi, I] = i{Kqi - X^Pi), 

\pi,Pj\ = il^^Fij, [qi, qj] = iX^Fij, 
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where Fij, pi, Xi are the generators of the Lorentz group and the operators of 
momentum components and coordinates, correspondingly, I is the "identity" 
operator, j, k,l = 0, 1, 2, 3. The new quantum constants fi and A have dimen- 
sionahty of mass and lenght correspondingly. The constant k is dimensionless 
in the natural system of units. 

By applying the algebra f[T^ to the description of color particles the con- 
dition K = can be imposed. Actually it is known the nonzero k leads to the 
CP— violation [31, 32], but strong interactions are invariant with respect to the 
P— , C— and T— transformations on the high level of precision. Moreover for 
color particles one can use the relation /lA = 1 [28]. In this case we obtain the 
reduction of GSYA to the special Snyder- Yang algebra (SSYA) with /xA = 1 
and K = for strong interaction color particles. Denoting /i as /ic and A as Ac 
we write the following commutation relations without the standard commuta- 
tion relations with Lorentz group generators, which are shown for the GSYA 
above (see eqs. f|T9|) ). 

[Pi, qj] = igiji, [Pt, I] = ifj'lqi, h, i] = -i^^lPh 

[qi,qj] = iXlFij, [Pt,Pj] = ifilFij. (20) 

We take into account difficulties arised when one try to prove the con- 
finement on the basis of the QCD first principles [33], and we simulate this 
phemomenon with the help of an assumed high symmetry of the nonpertur- 
bative OCD interaction beyond the Poincare symmetry. So we turn from the 
Poincare symmetry in the Minkowski spacetime to the inhomogeneous 0(2,6) 
symmetry in a phase space of a color particle [28]. 

Under these conditions the new Dirac type equation for a spinorial field ip 
has the following form: 

7^Pa^ = M^, 7^7^ + 7 V = 2g^^. (21) 

where 7"^ are the Clifford numbers for the spinorial 0(2,6) representation. 

One can take the product of Eq. f l2T|) with 7'^Pa + M and apply Eqs. f lT9|) . 
then the following equation for can be obtained 

{p'pi + a^q'q, + 2S,<,.S^^P,,+ 

+ 2a5°J)V = M^^, 5° = -CO, S'^ = -C'\ (22) 

where C°, C^^ are the combinations of the 7-matrices and the constants A, 
/i and K [28]. Eq. fl22|) contains the oscillator potential, which restricts a 



8 



motion of a color quark. Besides that we broke the inhomogeneous 0(2,6) 
symmetry with the help of the commutation relations fll9l) . In the special case 
fiX = 1, K = and the commutation relations fl20|) for SSYA we will obtain 
more simple expressions for the and but the form of the Eq. fl22|) will 
remain unchanged. Note that Eq. fl22|) can also be applied for a description 
of a confinement of boson particles such as diquarks and gluons with the same 
confinement parameter a. 

One can get an estimation of the a value using Eq. (!22|) . As it is seen, 
and p"^ entered into Eq. (l22!) can be considered as current and constituent 
quark masses squared, respectively. So Eq. fl22l) indicates that the convential 
relation = for a current quark should be transform to = + 
for a constituent quark, where = + is a constituent mass squared. 
To estimate the a value with the help of a constituent quark mass m and a 
current quark mass M values a ground state tpo in a meson has been considered 
neglecting the orbital angular momentum contribution LtpQ, as it had been 
done for DGLE. We can also estimate the a value with the help of the value 
of the confinement rising potential coefficient [3]. It is seen from Eq. fl22p 
that the coefficient of the oscillator confinement potential for a color particle 
is equal to o"^. Within the potential approach this coefficient is connected with 
the string tension astr typically as cr^ = where agtr varies from 0.19 

GeV^ to 0.21 Gel/2 [gj^ g^^^g ^ _ 0.3GeV and Xconf ~ O.QFm. Clearly it 
is assumed that the conceptions of the asymptotical oscillator potential and 
the constituent quark are applicable in a confinement domain. 

6. Conclusion 

The Dirac equation is successfully employed for the description of particles 
with the spin 1/2. If one take into account an additional Pauli term, then the 
Dirac equation describes in the SM framework the modification of a magnetic 
moments of leptons due to radiative corrections. However, someone has cast 
doubt on the relevance of DE for quarks in a confinement domain, i.e. in 
the nonperturbative domain of QCD. It is very likely that in this case we 
should modify essentially the Dirac equation or to pass on to a new equation 
of the Dirac type. The possible generalizations of the Dirac equation have 
been considered in this report. The Dirac equation in a quantum phase space 
is the most interesting from the standpoint of description of quark motion. 
This equation involves an oscillatory potential, so it can play a role of a model 
of color particle confinement. 

Ya.I. Azimov, S.V. Semenov 
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